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ABSTRACT 



During underwater (wet) welding the metal within the heat- 
affected-zone experiences a rapid fluctuation in temperature 
which results in a loss of ductility and a susceptibility to 
stress-corrosion cracking. The object of this study was to 
develop and verify a model which would predict the temperature 
distribution such that the metallurgical and mechanical pro- 
perties of the underwater weld could be determined. 

The equation of heat flow in two-dimensions with source 
terms was converted to its finite difference equivalent and 
solved by successive iterations with a digital computer. The 
solution incorporated all non-linear terms in the original 
equation and accounted for the variation of the thermal pro- 
perties of the metal with temperature. The interior boundary 
in the temperature field was found to have a definite shape 
whbh was a function of the important welding parameters. The 
shape function was correlated by experiment using methods 
established for normal surface welding. Heat losses were as- 
sumed to occur by boiling and radiation heat transfer. 

The predicted and measured temperature distributions were 
compared. It was found that the model predicted peak tempera- 
tures with reasonable accuracy but that comparison of cooling 
curves showed considerable disagreement. Since the predicted 
cooling curve presented a sharper slope than the measured value 
it was considered that the assumption of conventional boiling 
heat transfer processes was invalid. 
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CHAPTER I 
INTRODUCTION 

The object of this thesis was to develop a model of an 
underwater welding process such that the temperature history 
of the metal could be predicted. Particular attention was 
given to temperature changes within the heat affected zone. 

The model was constructed with consideration for a) the varia- 
tion with temperature of the thermal properties of the base 
metal, b) a finite heat source of predictable configuration, 
c) the heat losses to the liquid environment, and d) the non- 
linear components of the heat flow equation. Three major 
areas within the literature were investigated: a) under- 

water, b) normal surface welding, and c) boiling heat trans- 
fer . 

This study was conducted in recognition of the poten- 
tial benefits that could be realized if underwater welding 
were a viable process. Predicted trends in ocean exploration 
and exploitation indicate the use of large structures that 
could only be fabricated or repaired on site within the 
liquid environment. Even today, considerable savings could 
be realized if ship hulls, off-shore drilling platforms, buoys 
and salvage work were constructed or repaired without ex- 
plicit dependence on the requirement that the welding opera- 
tion be conducted in air. Unfortunately, underwater welds 
are not structurally sound due to reduced ductility, and 
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susceptibility to stress corrosion and hydrogen cracking. 

A reasonably accurate model, one that predicts peak 
temperatures and cooling rates, was seen as a necessity to 
examine the metallurgical changes occurring in the weld 
zone as function of welding parameter. 



A. UNDERWATER WELDING: STATE-OF-THE-ART 

Underwater welding, as a process, is best defined by 
considering its two subcategories: wet and dry welding. 

Dry underwater welding is a welding procedure in which the 
liquid environment is excluded from the vicinity of the work- 
piece by a structure such as a caisson. The structure is 
large enough to contain one or more workers plus the neces- 
sary welding equipment. The atmosphere within the chamber is 
usually at depth pressure and may, or may not, consist of 
life-supporting gases. This system's utility is restricted 
to symmetrical workpieces about which the chamber can be 
fitted and sealed. Hence, dry underwater welding is most 
used by oil companies laying undersea pipeline or repairing 
the submerged supporting structure on off-shore drilling 
rigs. Essentially, this procedure differs from normal air 
welding only by the pressure of the atmosphere surrounding 
the work. 

Wet welding occurs when the liquid environment is al- 
lowed to interact directly with the joining process: the 
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most notable result being a rapid quenching of the liquid 
metal. The utility of this process is not restricted by 
the geometry of the workpiece, but rather by the quality of 
the product. No wet weld is considered safe when used in 
load bearing structures. The rapid quench produces a weld that 
has 50% of the ductility of air welds, and is subject to 
extreme stress corrosion cracking and hydrogen embrittle- 
ment. The object of the thesis, then, was to construct, 
and verify, a model that might predict the temperature cy- 
cles in a metal as it is being wet welded. 

The necessity for such a model is twofold: a) considera- 

ble time and expense could be saved if current ocean struc- 
tures could be repaired or fabricated in situ , and b) future 
projections of man's involvement with the ocean require the 
use of structures that can only be fabricated, by reason of 
their size, within the ocean environment. Currently, the 
most economical means of joining metal is by welding, and it 
would be extremely expensive if each new process, each new 
structural material, with their multitude of variables, had 
to be investigated individually by actual experimentation. 
Unfortunately, the state of wet welding has generally stag- 
nated since its inception during WW II. 

From its development in the 1940's until today, the 
most widely used wet welding process is shielded metal-arc. 

This is the same method used in air welding and is commonly 



12 



referred to as "covered electrode" or "stick electrode" 
welding. Generally, wet welding has merely borrowed what- 
ever equipment was already developed for air welding. The 
only adaptation to the new environment being a coat of water- 
proof material for the electrode and an electrically-insulated 
holder. Only recently have the British begun to develop 
electrodes specifically designed for use under water. How- 
ever, the details of the program are considered proprietary 
information. 

The basic research conducted thus far has been an ex- 
amination of the mechanical properties of the flux and the 

23 

completed weld. Thus, Waugh and Eberlein investigated 
flux coverings to determine those most suitable for under- 

2 6 

water application.* Work at the Naval Engineering Station 
produced some factors of safety that must be applied when 

g 

wet welding. Hibshman and Jensen evaluated the properties 
of welds in low-carbon steel, under water and in air. The 
data revealed that ductility was markedly reduced. As far 
as it is known, however, no attempt has been made, either on 
a theoretical or experimental basis, to determine the tem- 
perature history of the metal during the wet welding process. 



B. NORMAL SURFACE WELDING 

A review of the literature for air welding, however, 
indicated that considerable effort has been devoted to 



* 



Numbers refer to listing in Bibliography. 
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predicting the temperature cycles during the welding process. 
In 1935, Rosenthal^ obtained an exact solution to a modi- 
fied heat flow equation. He simplified his model by neglec- 
ting non-linear terms, treating the conductivity as constant, 
and assuming a point source. However, the resulting equation 
for temperature distribution contained all the important weld- 
ing parameters: weld speed, heat input, plate thickness, and 

the thermal properties of the metal. Furthermore, the equa- 
tion was in quasi-steady state form so that the location of 
the point source was at the origin of a stationary coordinate 
system. Experimental verification by Rosenthal‘S and 
Schmorber indicated that the model was reasonably accurate 
at some distance from the source of heat, but in adjacent 
areas (within the heat affected zone, or haz) the degree of 
error was unacceptably large. The primary cause for the 
discrepancy was due to the nature of the heat source which 
in actuality has a very definite shape. Nevertheless, 
Rosenthal's initial attempt provided the impetus which spanned 
a succession of attempts to include all the important varia- 
bles in one well organized model. 

7 

Grosh, Hawkins, and Travant investigated the effect of 
considering the thermal properties as a function of tempera- 

3 

ture and obtained slightly better results. Barry modified 
Rosenthal's model by including possible heat losses to the 
environment but concluded that, at least for convection, 
the effect is negligible in comparison to the heat conducted 



14 



in the metal. 

The most important attempts to derive a satisfactory 

model concerned those which examined the shape of the heat 

27 

source. In 1952, Wells proposed an equation from which 
the length and width of the weld pool could be determined 
as a function of the heat input. Subsequent experiments 

indicated that the relationship was reasonably accurate. 

. . 4 

Christensen, Davies, and Gjermundsen examined the shape 

of the weld pool and its relation to a wide variety of 
welding conditions and associated variables. While their 
objective was not to formulate a model, they did present 
experimental data from which viable conclusions were drawn: 
a) pool width and length increase with heat input for a 
given weld speed and plate thickness, b) the pool is wider 
than predicted by point source theory, and c) the pool has 
a distinct "tear-drop" shape. These basic observations per- 
mitted the derivation of a very successful model of tempera- 
ture distribution during welding. 

13 20 

In 1967, Pavelic and Tanbakuchi constructed and 
verified a model which combined all the disparate elements 
necessary to predict temperature variations in the haz . 
Unlike Rosenthal, they considered the variation of conduc- 
tivity and diffusivity with temperature, a finite heat 
source, heat losses to the environment through convection 
and radiation, non-linearites in the heat flow equation, 



and preheating of the metal from arc radiation. Their ap- 
proach to the derivation of the model for a tungsten-inert- 
gas welding process formed the basis for this thesis. 



C. BOILING HEAT TRANSFER 

The literature on boiling heat transfer was investiga- 
ted selectively with respect to possible models relating 
to sub-cooled pool boiling on horizontal plates at a non- 
uniform temperature. 

The complexity of the boiling process has prevented 
researchers from developing a single formulation covering 

the wide range of temperatures over which boiling will oc- 

. 15 

cur. Figure 1-1, reported by Rohsenow, is considered to 

represent the regimes of boiling, each regime having a 

uniquely formulated model descriptive of heat transfer. 

15 

Rohsenow has correlated the data of Addoms for nu- 
cleate boiling on platinum wire in water at saturation tem- 
perature. Although the data was given for a specific ge- 
ometry, Rohsenow observed that the same results applied to 
horizontal surfaces. 

Correlations in the transition regime have not been 

25 

firmly developed. Jordan, reporting the data of Drew 
and Mueller, Farber and Scorah, and Braunleh, concluded that 
if the data were plotted on a log-log graph, they would lie 
along a straight line from the critical flux (or, burnout) 
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FIGURE 1-1 
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to the minimum film boiling heat flux. The slope of the 
line varied with the surface finish: the steepest slope oc- 

curring with a mirror smooth finish. 

25 

Chang, Botushansky, and Berenson developed correla- 
tions for film boiling heat flux from horizontal surfaces. 

25 

Jordan has reported that the Borishansky relation fits 
the data for pool boiling in water. 



D. SUMMARY OF PROCEDURE 

The development of a model for predicting temperature 
changes in the heat affected zone of an underwater weld 
required the construction of a theoretical model followed 
by experimental verification. The model was to encompass 
the latest techniques, and should be reasonably accurate 
in determining peak temperatures and cooling rates. 

The derivation of the model followed the approach taken 
13 20 

by Pavelic and Tanbachuchi in their work for surface GTA 
(Gas Tungsten Arc) welding. The two-dimensional equation of 
heat flow with external sources was reduced to quasi-steady 
state form and transformed to its finite difference equiva- 
lent. The external heat sources were examined and the proper 
relations chosen. Then, with the help of Pavelic's results, 
the equation describing the shape of the molten weld pool 
was derived. The resulting derivation included numerical 
constants which were determined through experimentation. 
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Finally, the computer program was written. The method for 
laying out the finite difference grid was that developed by 
Tanbachuchi. His solution was refined by having the shape 
of pool internally generated instead of being supplied as 
input . 

The experiment was conducted in two parts. First, the 
correlation constants in the pool configuration equation 
were determined by measuring the dimensions of molten weld 
pools while varying the welding parameters: Heat input 

and welding speed. Then, temperatures occurring during 
welding were determined by laying weld beads on submerged 
plates in which thermocouples had been installed. The re- 
sults of the final experiment were compared with the values 
predicted by the model . 
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CHAPTER II 

MATHEMATICAL MODEL AND COMPUTATIONAL TECHNIQUE 
A. PRELIMINARY OBSERVATIONS 

Welding is a process of joining two or more materials 
together by applying sufficient quantities of heat to a lo- 
cal area at their interface such that melting occurs and a 

union is achieved. The amount of heat required is on the 
4 

order of 10 BTU/HR and, for structural applications, the 
manner of generation is by electric arc. The method for 
directing the arc is by an electrode which may (consumable) 
or may not (non-consumable) contribute metal to the finished 
bond. The electrode concentrates the arc and creates a small 
pool of molten metal which is typically "tadpole" or "rain- 
drop" shaped when viewed in the plane perpendicular to the 
arc axis. As the electrode is drawn along the material a 
bead of weld (i.e., a chain of overlapping pools) is formed 
which represents the physical evidence of the interface bond. 

The largest percentage of the heat developed by the arc 
goes into raising the temperature of the material to its 
melting point. The remainder is accounted for by conduction, 
convection, and radiation, joule heating, and pre-heating of 
the material in the vicinity of the arc. Pre-heating will 
occur in all but submerged arc processes and is due to the 
divergence of the arc which takes the appearance of a cone 
with the electrode at its apex. The central part of the 
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arc, or the plasma column, is responsible for the creation 
of the molten pool. The remainder of the arc, while not 
generating sufficient heat to melt the material, will cause 
preheating. Hence, heat may flow into the material direc- 
tly from the pool and from the arc by pre-heating. 

While the above heat balance has been observed in nor- 
mal sruface welding, its complete relevance to the under- 
water welding process is less certain. Conduction and radia- 
tion have importance, but convective losses must be modified 
by boiling heat transfer phenomena. Joule heating losses 
are usually estimated to account for 1-2% of machine power 
and, hence, may be neglected. A search of the literature 
did not reveal any information on the dynamics of an arc 
surrounded by water so that the efeect of pre-heating could 
not be determined. For the experiment, it was decided that 
the effect, if it did exist, could be reduced to negligibi- 
lity by the proper choice of electrodes. Consequently, an 
electrode that could be placed directly on the metal sur- 
face without shorting out was selected. These "drag" type 
electrodes operated in a manner such that the metal core 
melted faster than the flux covering the flux coat, which 
came in direct contact with the metal, formed a hollow bar- 
rel and prevented arc diffusion. Consequently, it was assumed 
that pre-heating was negligible and all heat flowed ii\to 
the metal through the molten pool. Interaction with the 
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environment was assumed to occur by radiation and boiling 
heat transfer. 

An approximation to the real welding process was then 
seen to be a boundary-value problem in two-dimensions with 
negative source terms. The interior boundary in the field 
was viewed to be formed by the contour of the molten pool; 
the exterior boundaries were considered as the four edges 
of a rectangular surface. The pool contour temperature was 
known to be the melting temperature of the material with 
the edges at ambient. Radiation and boiling heat transfer 
were included as the negative sources. 

B. EQUATION OF HEAT FLOW 

The governing differential equations are presented in 
this section. In the sections to follow the equations are 
reduced to their finite difference equivalents, the boun- 
daries are defined, and the relations for the source terms 
are determined. 

The general equation of heat flow with source terms 

is, 

V- (K T) + W i = , 

where k, P , and c are the conductivity, density and specific 
heat of the material and may be a function of distance, time, 
and temperature. VT represents the source terms and has the 

3 

units, BTU/FT-HR. In partial differential notation for two 
dimensions, the equation becomes, 



The equation in this form is for a coordinate system 
fixed to the control surface defining the temperature field. 
Since the interior boundary (i.e., the pool contour) moves 
through the field with velocity, V, the velocity of the elec- 
trode, the temperature at a given point will change with 
time. To eliminate the time dependence, the coordinate sys- 
tem was transferred to the moving boundary with the follow- 
ing transformation equations: 



£ = x-Vt 

£ = t 



where £ and £ arc coordinates in the new system. (The de- 
tails of the calculation appear in Appendix 1.) With the 
transformation, the equation of heat flow became, 



2.1 



9 2 T , 9 2 T x 1 ak 

dx T By 7 + k 9T 




W. - V 9T . 
k a 9x 



The temperatures with respect to the new system could 
now be considered constant and were viewed to form isother- 
mal contours in the plane of the field with the pool contour 
forming the innermost. Hence, the pool boundary could be 
called the isotherm of melting. 



The equation was put in non-dimensional form by 
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defining normalizing relations. It was decided to make the 
exterior boundary conditions homogeneous and the conditions 
along the pool contour unity. Hence, the normalized tempera 
ture was defined, 



T* E 




2.2 



where is the melting temperature of the material (ap- 
proximately 270Q°F for mild steel) and T q the ambient 
temperature . 



The normalized dimensions were defined to be, 



x* = 
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y* = 



y_ 

x 

o 



when these definitions were substituted back into equation 

2.1, it was found that x_ should be defined as 

o 

x q = & /V 

where d refers to a specific value of the thermal diffu- 

sivity. From a stability analysis of equation 2.1, recom- 
2 

mended by Ames, it was determined that 



where aMin represents the minimum thermal diffusivity for 
the material. Finally, the normalized dimensions became 



x* = 



XV 

2 a Min 



T* = 



YV 



2a 



Min 



2.3 
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by substitution of equations 2.2 and 2.3 into equation 2.1, 
the heat flow relation was obtained in non-dimensional form, 

2 c( 

Min 3T 
a 3x 

2.4 

where the "Star" superscript has been dropped. (See Appendix 
4 for details of the derivation.) 



3 2 T 3 2 T ± 1 dK 

Jx* k dT 



(91) + (9T) 

^ 3x' ^ 3y ' 



. 4 aMi n j 
KV 2 (T m -T q ) 1 



B. FINITE DIFFERENCE SOLUTION 



Equation 2.4 was then converted to its finite difference 
equivalent so as to make it amenable to computer solution. 
Following the procedure recommended by Greenspaif for elliptic 
non-linear partial differential equations. A finite differ- 
ence grid node (Figure 2-9) was defined and the bounding tem- 
peratures were approximated by a Taylor expansion about T q : 



T. = T (±) s i 3T ° + s i 

1 ° ~5x 2~ 



2 



2m -i 



3^T 

die 



i = 1,2, 3, 4 



Now, in the limit as the S^ (i.e., the grid spacing) approach 
zero the variation in temperature between the bounding tem- 
peratures along a given axis will represent the actual rate 
of change of temperature at the node. For unequal grid spac- 
ing, these changes were defined as: 



3T 

tt — = a _T A + a.T, + a 0 T 0 

3x 00 11 22 



.3 t 
3y 



(2 ' 5 a,b> 



- = a T + a. T + a T 
0 o 1 3 2 4 
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The weighing coefficients were found by comparison with 
Taylor's expansion through combination of like terms. The 
second partials were similarly defined 



8 2 T 



3x 
8 2 T 



■z^- = a 3 T 0 + a 4 T x + a c T 



■ay 



5 2 



<2 - 6 a , b > 



■ 2 — = ^ 3 T 0 + ^ 4 T 3 + ^ 5 T 4 



and the coefficients found from Taylor's expansion. Substi- 
tution of equations 2.5 and 2.6 into 2.4 resulted in the finite 
difference equivalent for the equation of heat flow: 






T. =< 



2a + S 2 i 



aS ^ ( S^+S 2 ) 



T 

+ 2 



2a - S. 



aS 2 ( S ^+S 2 ) 



+ T 3 



s 3 (s 3 +s 4 ) 



+ T 4 



S 4 (S 3 +S 4 ) 



4a 



Min 



+ KV 2 (T m -T ) [W i J 
M o 



1 dK 
k dT 



S -S 

[( 1 



S 1 S 2 } T ° + ( S 3 (S 3 +S 4 ) 1 +( S 2 (S 1+ S 2 ) } 2 



2 



i /s S) ° +( 

b 3 b 4 



-) T, + (: 



S 3 (S 3 +S 4 ) 7 "3 ' s S 4 (S 3 +S 4 ) 



-) T 



>/' 



2.7 



where , 



B = ( 



S 1 S 2 



S 3 S 4 



S 2' S 1 

+ ( = Z ) 

^a(S 1 S 2 ) J 



The details of the derivation are presented in Appendix 
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Equation 2.7 operated on a grid (i.e., an assemblage of 
nodes) placed over the temperature field by assigning to each 
node point a value consisting of the weighted "average" of 
the four bounding temperatures. Since the exterior points 
in the field were at a normalized value of zero it was de- 
cided that calculations should begin at the pool contour 
(T - 1) and move outward. A series of iterations over the 
grid were performed until the desired accuracy was obtained. 



C. THE BOUNDARY CONDITIONS 

This section deals with the derivation of the equation 

describing the shape of the pool contour (i.e., the interior 

boundary in the temperature field) and a determination of the 

location of the exterior boundaries. The treatment of the 

13 

problem follows the approach taken by Pavelic in his study 
of the tungsten inert gas process for normal air welding. 
Pavelic' s method for finding the shape of the pool can be 
summarized as follows: 

1. Rosenthal's two-dimensional equation with a point 
source of heat is analyzed to obtain the relations 
predicting the maximum dimensions of the pool: 
length and width. 

2. These theoretical relations are correlated with 



experimentally-obtained data. 
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3. The contours of actual pools are measured. and a 
curve fitted equation is obtained. This equation 
is a function of the maximum pool dimensions. 



Hence, for any given set of welding variables, Pavelic could 
predict the contour of the pool. The details of Pavelic' s 
analysis have been summarized in Appendix 3. 



It was found that the maximum theoretical dimensions 
of pool are 



\ 



Y = r 
max m 



1 - 



V r m> 

K 1 (r m> 



i 






For the maximum half width, and 



X = 2r K o (r m^ 

max m K (r ) 
l m 



for the maximum length. The argument of the Bessel functions 
K q and , is the radial distance from the origin of the arbi 
trary coordinate system to the pool contour and was found to 
be a function of the lumped parameter, p: 

3.415 El 
p 2 tt a L H 

m. 

The relationship of x , Y , and r to the pool is pre- 

* max max m * c 

sented in Figure 2.1. When plotted on a log-log graph, the 
pool width (= 2y^ ax ) vs p and x max vs P appear as a straight 
line and, hence, experimental data could be correlated by 
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the relations: 



2y 

J max 

x 

max 




with a, b, c, and d becoming the correlation constants. 

Experiments were conducted therefore to find values 
for the constants. Pools were made by welding beads on 
thin rectangular plates. At a predetermined point in the 
welding process the molten metal was ejected from the pool 
by a jet of high pressure air leaving a crater whose edges 
formed the desired contour. Only those craters nearest 
the ideal two-dimensional configuration were retained for 
measurement. An ideal 2-D crater has been defined in 
Figure 2.2. The maximum length and width were measured 
with a Bausch and Lomb microscope containing a movable 
micrometer table. A computer program was then written to 
fit a curve through the data points from which values for 
the correlation constants were obtained. The program with 
sample input is presented in Appendix 7. 

Two major problems were recognized early in the attempt 
to produce craters for measurement: a) a search of the litera- 

ture did not reveal any similar experiments conducted with 
stick electrodes, and b) it was very difficult to time the 
release of air to eject the metal during the under water 
welding process. Without the certainty that the metal was 
being ejected at the proper instant it was considered 
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possible that the craters might be distorted and measure- 
ments valueless. The lack of prior experiments conducted 
in air denied a basis for comparison to determine if the 
underwater data were reasonable. Consequently, it was de- 
cided to conduct an initial series of runs in air where 
crater production could be closely controlled, merely to 
serve as a crude reference for the underwater results. 

Plates of two different thicknesses, 1/8" and 1/4", 
were used and the data plotted as shown in Figures 2.3 and 
2.4. The results for the 1/8" plate appeared to be in 
good agreement with the predicted relation, but those for 
the 1/4" plates were not. The lack of alignment was attri- 
buted to the failure to obtain two-dimensional pools in 
the larger specimens. The available power was insufficient 
to create the required configuration, and, as a result, the 
crater formed was only a hollow depression on the surface of 
the specimen. Significantly, an increase in power (ExI) 
would have resulted in a larger value for the lumped parame- 
ter, p, which would have shifted the data points to the 
right and brought them more in line with the predicted 
relationship. (It was noted, however, that a line drawn 
through the 1/4" data points was approximately parallel to 
the theoretical line. This relationship suggested that the 
penetration of the arc could also be correlated with the 
lumped parameter, p, so that the equation for a 3-D pool 
might eventually be derived. No attempt was made to do so, 
however.) From the plot for air welding it was observed 
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that a) 2-D pools would not be obtained in water for 1/4" 
thick specimens and b) the width of the underwater pool 
should be larger than predicted while the length should 
be smaller. 

Pools were then created under water, measurements taken 
and the data plotted as shown in Figures 2.5 and 2.6. The 
location of the data points with respect to the theoretical 
line was noted and compared with air welded results. Since 
the location was as expected it was assumed that the measure- 
ments were reasonable. Curves were fitted through the data 
points. These curves were not found to be parallel either 
to the predicted or the air welded lines, and the reason for 
the divergence (considered more fully in Chapter IV) could 
not be adequately explained due to the small size of the 
data sample and data scatter. Nevertheless, the correlation 
constants obtained from the fitted curve were used. It was 
consequently assumed that the length and width of an under- 
water pool were related to the lumped parameter, p, ac- 
cording to the equations: 



x* = 4 . 4 2 p 
M r 



1.01 



W* = 3.16 p 
M 



0.67 



where x£ and Wj* are dimensionless quantities and 
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with the maximum pool dimension relations defined it was now 
possible to derive an equation that would predict the com- 
plete contour of the pool. Pavelic suggested curve fitting 
an equation of the form: 

Y 2 = Bx a (x m - x) b 2.8 

to which the boundary conditions were applied, 
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where a and b are the correlating constants. The details of 
the derivation have been included in Appendix 3. 

Values for the constants were obtained by measuring 
the actual pool contours and fitting the resulting measure- 
ments with equation 2.8. A computer program, employing an 
M.I.T. non-linear least squares curve fitting routine, was 
written and has been presented in Appendix 7 . It was hoped 
that a large number of pools would be available for measure- 
ment in order that the calculated correlating constants 
would be generally applicable. While a large sample size 
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was obtained for air welded pools, only one under water pool 
of acceptable quality was produced. Due to the accelerated 
solidification of the molten metal in water and crater walls 
took on a jagged appearance which made measurements diffi- 
cult although the maximum dimensions were generally easily 
determined. Consequently, it was decided to base the cal- 
culation of the constants on data taken from the air welded 
pools and then compare the predicted contour with the mea- 
surements from the single under water pool. The result of 
this approach is presented in Figure 2.7. The predicted 
contour gave a wider head and narrower tail than was mea- 
sured but the overall shape was considered satisfactory. 
Nevertheless, it was realized that the contour equation was 
based on a sample size of one (i.e., the single under water 
pool) . However, it seemed reasonable to assume that the 
general shape of a pool would be independent of the atmos- 
phere in which it was created sirce the dimensions were most 
strongly influenced by heat input and welding speed. With 
this assumption, the correlation constants 

a = 0.83 
b = 2 . i.6 

were used to predict the shape of the under water pool con- 
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D. DERIVATION OF THE HEAT LOSS TERMS, W. 
' ! 

For this model heat losses were assumed to occur 
through boiling and radiation. Although convective heat 
flux could be considered, it was estimated that these losses 
would account for only 3% of the total heat input to the 
specimen (15,000 to 20,000 BTU/Hr) and were therefore neg- 
lected. Consequently, heat losses were considered accord- 
ing to three distinct regimes as depicted in Figure 2.8. 

Nucleate boiling is characterized by the termination 
of bubbles at the heated surface which grow in size with 
increasing temperature. The agitation of the liquid near 
the surface caused by rubble activity — formation, collapse, 
or ascent — maintains a circulation of relatively cooler 
liquid near the surface which can result in heat fluxes on 
the order of 10 6 BTU/FT 2 -HR. 

When the surface temperature reaches a value 50-90°F 
above saturation the number of bubbles forming has increased 
to the point where they begin to overlap and create a solid 
film of vapor. The vapor film is not immediately stable and 
may collapse back to the heated surface where the process 
repeats itself. In this transition region the cooler liquid 
is brought into contact with the surface only intermittently 
and heat flux rates are not as high as they would be during 
nucleation. As temperature increases the film grows more 
stable, insulating the surface and reducing heat transfer. 
Finally, the film becomes fully stable and the heat flux 
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reaches its minimum value. At higher temperatures heat loss 
through radiation becomes the predominant term. 

The complexity of the entire process required that each 
regime be treated separately in determining the appropriate 
relation for calculating heat flux. First, however, it was 
necessary to define the temperatures which bounded each re- 
gime. Since nucleation begins near saturation temperature, 
the lower bound was chosen to be 212 °F. Burnout, or the end 
of nucleation, will occur at a temperature which depends 
upon the condition of the heated surface, the amount of dis- 
solved gases in the liquid, and the ambient pressure. If 
the solubility of air at 14.7 psia is 6m SL/SL, then, from 
data reported by Rohsenow^ burnout will occur neat 260 °F. 
However, a mildly oxidized surface increases the tempera- 
ture, by an unknown amount. Consequently, 300 °F (as the 
nearest round number) was defined as the temperature at 
which nucleation terminates and the transition regime begins 

For the minimum heat flux, the data of Borishansky and Beren 

25 

son as reported by Jordan indicated that transition ends 
near 500°F. 

Once the boundary temperatures for each boiling regime 
were established it was necessary to determine the suitable 
heat flux relations. Rohsenow's^ 5 correlation of Addoms ' 
data suggested that heat flux in the nucleate boiling regime 
varied according to the cube of the temperature difference, 
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The properties of the liquid and vapor were calculated at 
saturation and the proper relation was found to be 

q N = [1.54 2 (T w -T sat )] 3 , BTU/ft 2 -hr. 

The details of the calculation have been presented in Ap- 

25 

pendix 2. As indicated by Jordan , the data for film boil- 
ing on a horizontal surface in water is predicted withrea- 
sonable accuracy by Borishansky ' s relation. In this regime 
the heat flux was found to be 

q F = 32.4 (T w - T sat ) 0 * 75 , BTU/hr-ft 2 

where calculations have been summarized in Appendix 2. 

Heat flux during transition boiling has not been well 
defined. However, data reported in the literature indicates 
that measures fluxes tend to fall along a straight line 
plotted with log-log coordinates. This relationship sug- 
gested the derivation of an equation with the burnout and 
minimum film boiling heat fluxes as boundary conditions. As 
summarized in Appendix 2, the results of the derivation were 
found to be 

q T = 2265.11 

where 

300° - T tt - 500 °F 

w 



288 



(T W t sat^ 



5.86 
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The equation for radiation heat transfer was determined 
by assuming the surface of specimen exchanges heat with the 
surrounding water. Hence, 

q R = 0.1317 x 10" 8 (T w 4 - T sa ^) , BTU/hr-ft 2 
E. COMPUTATIONAL TECHNIQUE 

The computer program written to calculate the temperature 
distribution consisted of the four main operations: 

1. The shape of the pool contour is generated. 

2. The grid is layed out to form a computational block. 

3. The grid nodes are initialized using Rosenthal's 
temperature distribution equation. 

4. Successive iterations are carried out over the 
grid until the desired accuracy is obtained. 

A listing of the program is presented in Appendix 7, and 
instructions on its use in Appendix 6. 

Pool Contour Generation 

A value for the lumped parameter, p, is calculated from 
the input data and the maximum dimensions of the pool are 
determined according to the correlation derived in Chapter II. 
Then, using the contour equation, which is a function of the 
maximum dimensions, the computer generates a series of points 
which define the shape of the interior boundary in the tem- 
perature field. If projected on the x-axis (i.e., the 
longitudinal center line) of the pool, the points would be 
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spaced 5/100 of an inch apart. This spacing allows from 10 
(0.5 inch pool length) to 20 points in defining the shape of 
the boundary. 

Grid Layout 

Except for minor alterations, the method of arranging 

. . . 20 
the grid is identical to approach developed by Tanbakuchi 

and is sized approximately 100 lines to the inch in a recta- 
linear array. The grid is configured such that the specified 
thermocouple location appears at the intersection of a row 
and a column. A schematic of the grid forming a computa- 
tional block is presented in Figure 2.9. Columns are ar- 
ranged parallel to the center line (i.e., column number 1) 
of the pool; rows perpendicular. The intersection of column 
1 and row 1 mark the origin of the coordinate system. In 
the output, rows forward, or in the direction of welding, 
are identified as positive quantities. 

Initializing the Grid Nodes 

Nodes within and on the pool boundary are given values 

of one (T* = 1) while the exterior boundary nodes are set to 

20 

zero. As suggested by Tanbakuchi, the intervening nodes are 
initialized with Rosenthal's temperature distribution equa- 
tion. Initializing the nodes considerably reduces the number 
of iterations required to define the temperature field. 

Calculation of Temperatures 



A Gauss-Seidel type iteration, which makes use of the 
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latest available temperatures, is performed until the speci- 
fied accuracy is reached. The program begins each iteration 
cycle at the interior boundary to avoid excessive computation 
in that part of the field initially at zero normalized- 
temperature. Calculation of the value for the heat loss term 
is accomplished by a search routine which locates the node 
temperature on the boiling curve and computes the appropriate 
quantity . 



48 

CHAPTER III 

UNDER WATER EXPERIMENT 

The experimental verification portion of the thesis 
consisted of three steps: Production and measurement of 

molten pools made in air, production and measurement of 
molten pools made under water, and measurement of tempera- 
ture changes in a specimen being welded under water. 

A. EQUIPMENT 

The same equipment shown schematically in Fig. 3-1, 
was used for all three steps mentioned above and consisted 
of the following: 

1. Water Tank 

2. Gravity Welding Stand 

3. Power Source 

4 . Blow-Out System 

5. Temperature Recorders 

6. Current, Voltage, and Time Recorders. 

Water Tank. The welding experiments were conducted in a steel 
tank measuring four feet long by three feet high by two feet 
deep. The exterior of the tank was sheathed in wood and rub- 
ber matting to prevent electrical shocks. A platform, con- 
taining clamps for securing the specimen, was fitted to the 
inside of the tank. The platform was located one foot from 
the top of the tank so that the electrode, but not the elec- 
trode holder, would be submerged in the water. The tank was 
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filled by a line supplying ordinary tap water (Boston, 

Mass.) and drained by gravity discharge. The tank was grounded 
to the power source. 

Gravity Welding Stand . A Japanese gravity welder was modi- 
fied and fitted to the top of the tank. The device consisted 
essentially of a common stick electrode holder fastened to a 
carriage which could be moved along an inclined rod. A guide 
was fitted to the rod so that the electrode would draw a 
straight bead on the specimen. The holder was adjusted to 
the carriage until the electrode maintained a constant angle 
of fifty-one degrees relative to the surface of the specimen 
in the plane of travel and ninety degrees in the plane per- 
pendicular to direction of electrode travel. The tripping 
mechanism on the carriage, which rotated the electrode up 
and away from the specimen at the end of travel, was adjusted 
so that the arc would be broken after approximately six in- 
ches of weld bead was drawn. 

Power Source . A Lincoln arc welder capable of supplying up 
to 40 volts and 350 amperes direct current was used as the 
power source. The device was set to operate on reverse po- 
larity. This arrangement was selected as being similar to 
power sources for actual under water welding. 

Pool Blowout System . A system for ejecting the molten metal 
from the specimen to form a weld crater was constructed and 
consisted of a nozzle, foot switch, and nose connected to a 
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400 psi air line. The nozzle was mounted on a swivel-necked 
holder so that adjustments could be made in directing the 
flow of air into the molten pool. The holder and nozzle 
combination were clamped to the welding platform inside the 
tank. Several preliminary experiments were conducted to de- 
termine the direction and angle of the jet. After consider- 
ing several possible locations, it was decided that having 
the nozzle discharge along the axis of the bead but opposite 
to the direction of welding would give satisfactory results. 
The nozzle angle was adjusted so that the air would jet into 
the head of the pool, directing the molten metal back into 
the tank, away from the operator stationed behind the gravity 
welding stand. Once the proper relationships were determined, 
the nozzle was securely fastened in place for the duration 
of the experiment. 

Temperature Recorders. Three clevite brush, Model 220, re- 
cording voltmeters were used to measure the output signal 
from the thermocouples. Each instrument contained two chan- 
nels so that a total of six thermocouples could be recorded 
simultaneously. The event markers were connected electri- 
cally and repressing the key on any instrument would cause 
a simultaneous voltage signal to be recorded on all six chan- 
nels. The instruments were calibrated at 200°, 700°, and 
1000°F by inserting a thermocouple, of the same type and 
dimensions used in the experiment, into a controlled bath. 

The chart speed was adjusted to 25MM/sec and the gain to 
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lMV/Division (or, 50 MV total scale). 

Common mercury-in-glass thermometers were selected to 
determine water and air temperatures. One thermometer was 
taped near the terminals of the recorders to measure the 
cold junction temperature; the other was suspended in the 
water to the depth of the specimen (approximately six inches 
below the surface at the far end of the tank) . 

Current, Voltage and Time Recorders . A General Electric 
D-C (0.400 amp.). Model 10H, recording ammeter was used to 
measure the current flowing in the welding circuit. An 
Esterline Angus (0.100V) recording voltmeter was selected 
to measure the operating voltage during welding. A timing 
device consisting of a General Time Co., Model 320-1, elec- 
tric clock connected to a circuit breaker was used to mea- 
sure the amount of time required to lay a given length of 
weld bead. When the arc was initiated the breaker would 
close and start the clock; breaking the arc would cause the 
clock to stop. Dividing the length of the bead by the elapsed 
time thus recorded served to indicate the welding speed. 

All three instruments had been calibrated by the in- 
strument lab at the Boston Naval Shipyard and were assumed 
to be recording properly. 

B. MATERIALS 



The criteria used in selecting materials for the 
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experiment was based upon their being readily available at 
low cost and having wide usage in the industry. Three main 
types of materials were used: Steel plate, welding rod, 

and thermocouple wire. 

Steel Plate . All specimens were cut from cold rolled mild 
steel (A.I.S.I. 1010) sheets. Three different thicknesses 
were selected: 1/16", 1/8", and 1/4". For the pool measure- 

ment portion of the experiments the sheets were cut into 4" 
by 6" sections; the specimens in which temperatures would be 
measured were sized 6" by 12". Both sides of all specimens 
were cleaned by grit blasting to provide a uniform surface. 

Welding Rod . E7018-A1, 3/16" Dia., stick electrodes were 
found to be suitable for the experiment. This metal core of 
rod tended to melt faster than the flux covering disintegrated 
permitting the electrode to be "dragged" across the specimen 
during the under water operation. This feature was critical 
since the arc would not be observed while under water. In 
essence, the flux covering formed a hollow cylinder which 
shielded the arc from the water, and maintained the proper 
length necessary to sustain the arc. Since the duration of 
submergence in water was short (20 to 30 sec!), the electrodes 
were not coated with any water-proofing materials. 

Thermocouples . Chromel-Alumel Type K, .020" Dia., thermo- 



couples were selected to measure temperature changes occurring 
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during welding. The selection of the size required a com- 
promise between the ability to withstand high temperatures 
(large diameters) , and fast response (small diameters) . 
Furthermore, the thermocouple junction should form a point 
(ideally) in the specimen to avoid thermal gradients. 
Consequently, the diameter chosen was the smallest considered 
practical for the experiment. 



C . EXPERIMENTAL PROCEDURE 

The experiment consisted of three parts: a) production 

of pool craters, b) temperature measurement, and c) reduction 
of data. 

Production of Pool Craters . Molten pools were first produced 
by welding in air, then under water. The essential steps for 
this part of the experiment were the same for both processes: 

1. Specimen aligned on work table. 

la. If an under water pool was desired the tank was then 
flooded until the specimen was covered by six in- 
ches of water. 

2. The power source was turned on and given a five 
minute warm up period. 

3. An electrode was selected and inserted into the 
holder . 

4. The holder and carriage on the welding stand were 
manipulated so that the tip of electrode was slowly 
brought to the point where the arc would start. 
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5. The carriage was drawn manually down the incline 
of the welding stand, laying a bead on the surface 
of the specimen. 

6. The blow-out system was activated by depressing 
the foot switch at the same time the carriage 
reached the tripping point. 

7. The voltage, current, time, ambient temperature, 
and bead length were recorded for the run. 

8. The dimensions of the pool crater - maximum width 
length, and contour - were measured under a Bausch 
and Lomb microscope equipped with a movable, micro- 
meter table. For the experiment, the ranges of the 
variables involved were: 



Current 
Voltage 
Welding Speed 
Ambient Temp 



100 - 300 amps 
20 - 30 volts 

50 - 125 Ft/Hr 
70 - 73°F (air) 

40 - 43°F (water) 



Current was varied by adjusting the rheostat on the power 
source in steps of 25 amps., while the voltage was allowed 
to seek its equilibrium value. A particular welding speed 
could not be chosen prior to a run since manual control was 
required. However, an acceptable range of speeds were pro- 
duced by consciously moving the carriage along the guide 
rod. A mechanical system for better control of welding 
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speeds was attempted but the problem of starting the arc 
defeated the project. 

The problems connected with producing two-dimensional 
weld pools by stick electrodes were numerous and it is esti- 
mated that only 20% of the runs resulted in acceptable data. 
The welding stand did not allow the freedom of movement 
that is normally enjoyed in "Striking" an arc. The low volt- 
ages required so small a gap that the electrode was often 
grounded before the arc could start. In addition, a precise 
relationship between the heat input and welding speed neces- 
sary to produce two-dimensional pools was observed. Since 
the speed could not be predetermined, there was no assur- 
ance that an acceptable pool would result from a run. In 
practice, a number of runs at various speeds were conducted 
at a given heat input until the proper pool was obtained. 

Temperature Measurement . All runs were conducted under water 
and the essential steps involved were: 

1. Preparation of the specimen. 

2. Aligning the specimen with the electrode. 

3. Flooding the tank until the level of water was six 
inches above the surface of the specimen. 

4. Turning on the power source and recording equipment. 

5. Welding a bead on the surface of the specimen. 

6. Recording temperature, time, voltage, current 
and bead length for the run. 
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The preparation of the specimen involved drilling holes 
at specified locations and mounting the thermocouples. 

Twelve .040" Dia. holes were drilled in each plate in a man- 
ner that assured that the axis of the hole was perpendicular 
to the surface of the plate. (See Fig. 3.2 for relative lo- 
cations.) The centers of the holes formed a locus of points 
five degrees from the center line of the plate. An addi- 
tional hole (#12) was drilled near the edge of the specimen 
to record boundary temperatures. (During the experiment, 
however, it was observed that placing a thermocouple at this 
location -- about six inches from the center line -- was un- 
necessary since the temperature variation one-half inch 
from the center line was negligible!) . The distance from 
the center line to the center of the holes varied from 0.125" 
for thermocouple locations #1 and #9 to 0.440" for location 
#8; the center-to-center spacing along the center line of 
the specimen was approximately one-half inch. This arrange- 
ment permitted temperature measurement within the heat 
affected zone but allowed sufficient spacing for mounting 
thermocouples . 

The thermocouples required special preparation to pre- 
vent their shorting out in the water. As shown in Fig. 3-3, 
water-proofing required sealing all joints with expoxy 
cement and covering the leads with a plastic sleeve shrunk 
to form a tight fit. (In their commercial form, the thermo- 
couples had already been covered with asbestos and silicon 
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FIGURE 3-2 
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impregnated fiberglass.) A ceramic bead, containing separate 



holes for each thermocouple, was installed to prevent the 
leads from making a junction outside the plate. The exposed 
leads were cleaned with emery cloth and inserted through one 
of the holes in the specimen. A pair of pliers were then 
used to draw the leads through the hole until the ceramic 
bead was nearly flush with the bottom of the plate. The 
edges of the hole on the top surface of the specimen were 
given a light peening with a small prick punch, and the leads 
extending above the specimen were then cut off and ground 
flush with the surface. The purpose of the peening was to 
give the thermocouples a preferential junction at the sur- 
face of the specimen, and avoid, as much as possible, thermal 
gradients in the specimen resulting from pools that were not 
exactly two-dimensional. In this way, six thermocouples 
were installed in each specimen, leaving six open holes. 

Twelve specimens were prepared such that no specimen dupli- 
cated another with respect to thermocouple locations. 

All specimens, with one exception, were 1/8" thick: 
the 1/16" plates burned through too easily and the penetration 
required for two-dimensional pools was not obtained with 
the 1/4" plates. One 1/4" specimen was prepared, however, 
in order to determine the response of the model to a situa- 
tion where the shape of the heat source was not the ideal 
configuration . 

A prepared specimen, containing six thermocouples, 



61 

was then placed on the table in the tank, aligned .with the 
electrode, and clamped firmly in position. (It was first 



hoped that the electrode could be made to follow exactly the 
center line of the plate, thereby simplifying measurements 
taken to locate the thermocouples. Actual runs revealed, 
however, that although the bead would run straight it could 
not be made to follow any predetermined path with sufficient 
accuracy. Consequently, the positions of the thermocouples 
with respect to the longitudinal axis of the weld bead were 
re-calculated following the run.) The tank was filled with 
water until the specimens were submerged to a depth of six 
inches. The thermocouples were connected to the recorders 
and the power source turned on. The electrode was guided 
into position in preparation for starting the arc. At the 
instant of arc initiation the laboratory assistant depressed 
the event key on a recorder which placed marks on all six 
output tapes simultaneously. The carriage containing the 
electrode was then either allowed to move freely or was 
drawn by hand with a slight steady pressure down the incline 
of the welding stand until the tripping point was reached. 
The ambient temperatures (air and water), run time, current, 
voltage, and bead length were recorded. 



D. DATA REDUCTION 

The maximum dimensions of the weld pool were correlated 
to the lumped parameter, p, by means of a least squares curve 
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fitting routine available as an option from M.I.T.'s computer 
library. The details of the program have been presented in 
Appendix 7. The correlating parameters for the equation 
determining the contour of the weld pool were found by a 
non-linear least squares curve fitting routine. (Also see 
Appendix 7.) 

The Y-displacement of a thermocouple was measured along 
a line running through the center of the thermocouple hole 
and perpendicular to the center line of the weld bead. The 
X-displacement was measured along the center line from the 
point of arc initiation to the intersection of the perpendi- 
cular. A computer program was written to reduce the thermo- 
couple output. The program scaled the distances on the out- 
put tape by the ratio of the chart to the welding speed, 
and ensured that the measured temperatures were in the same 
coordinate system as the model. The program, with sample 
input data, appears in Appendix 7. 
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CHAPTER IV 

RESULTS AND CONCLUSIONS 

A. TEST RESULTS 

The experimental portion of the thesis consisted of 
two parts: a) determining the correlation constants used 

to define the shape of the heat source, and b) measuring 
the temperature distribution in the specimen during weld- 
ing. The results of the pool contour study have already 
been presented and described in Chapter II; it need only 
be added here that the results were in good agreement with 
theoretical predictions. 

The results of comparing the temperature distribution 
as predicted by the model with that measured by experiment 
is represented in Figures 4-1 and 4-2. The coordinates 
have been presented in dimensionless form. The ordinates 
may be transformed to real temperatures by the relation 

T = 2677T* + 503 , °R 

The abcissa for Figure 4-1 can be converted to inches 
with the relation 

x = x*/40.59 , in. 

and for Figure 4-2 

x = x*/41.02 , in. 

where the starred terms are the dimensionless values. 

In the figures, Y is the distance from the weld center 
3 m 



FIG 4-1: CALCULATED AND MEASURED TEMPERATURE CURVES 
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FIG 4-2: CALCULATED AND MEASURED TEMPERATURE CURVES 



3tfniVd3dH3i 9S31NOISN3HIC3 
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DIMENSIONLESS DISTANCE 
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line to the point of maximum pool width; (Figure 4-1) and 
(Figure 4-2) are the locations in the specimen at which the 
temperatures were measured. 

Figure 4-3 has been included so that a comparison might 
be made between the predicted temperature distribution for 
underwater and air welding: the curve labeled "without boil- 

ing" represents air welding. 

Figure 4-4 presents the predicted temperature distribu- 
tion at various distances from the center line. 

Figure 4-5 represents a comparison between measured 
and calculated cooling curve slopes as a function of distance 
from the maximum pool width. The slope of a given curve was 
measured at the average temperature and was defined to be 
the change in dimensionless temperature with respect to the 
change in dimensionless distance. For example, the average 
measured temperature for Run 104 (Figure 4-1) was found to 
be 0.4 and the slope at that point was -0.025. 



B. ANALYSIS OF DATA 

Prediction of peak temperatures within the haz was found 
to be reasonably accurate. The average error was 13% of the 
measured temperature with the largest variation occurring at 
locations nearest the molten pool. 

The cooling portions of the predicted and measured curves 
were not found to be in good agreement. As a quantitative 



FIG 4-3: TEMPERATURE COMPARISON CURVES 
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FIG 4-4: THEORETICAL TEMPERATURE DISTRIBUTION 
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DIMENSIONLESS DISTANCE 



FIG 4-5: COMPARISON OF COOLING CURVE SLOPES 
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indication slopes were measured at the average temperatures 
of the distribution curves and plotted as in Figure 4-5. 

Near the molten pool the predicted slopes were in error by 
as much as 200% and converged with the measured values only 
at the edges of the haz. 

Consequently, an evaluation of the effect of boiling 
heat transfer was made by neglecting such losses (i.e., as 
if the specimen were welded in air) with results presented 
in Figure 4-3. Such an analysis indicated that conventional 
boiling models were inappropriate when applied to the under- 
water welding process. However, nothing in the data suggested 
what mechanism might be controlling heat losses. 

It can be speculated, nevertheless, that the gases gen- 
erated from the oxidizing flux and possible dissociation of 
water molecules could disrupt the normal process of boiling 

i 

heat transfer. These gases would collect at the tip of the 
electrode to form a protective envelope over a portion of 
the specimen's top surface. Within the envelope the predomi- 
nant form of heat transfer would be radiation while exterior 
areas could support some boiling losses. The bottom surface 
(opposite to the side being welded) would not see the shield- 
ing gases and, consequently, heat losses would occur through 
boiling. This mixed form of heat losses would explain the 
slower rate of cooling observed during the experiment. 

Unfortunately, the process just described has not been 
described in the literature and the development of the 
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appropriate model was beyond the scope of this research 
study . 



C. CONCLUSIONS 

The object of this thesis was to develop a model that 
would predict the temperature distribution in thin plates 
welded underwater. 

It was found that the model was reasonably accurate 
in predicting peak temperatures within the heat-affected- 
zone, but that the degree of error in predicting cooling 
rates was unacceptably large. The initial assumption that 
heat losses can occur by conventional mechanisms of boiling 
heat transfer was, therefore, considered invalid. Rather, 
the proper model for heat losses must consider the effect 
of the gases generated during the welding operation. 

Given the above conclusions, it is recommended that 

1. Heat transfer fluxes be correlated for the 
case involving full boiling with radiation on 
one surface while the opposite surface supports 
only radiation, and 

2. The study should consist of vertical as well as 



horizontal surfaces. 
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APPENDIX 1 

TRANSFORMATION TO QUASI-STEADY STATE 



K 



9 2 T 
dx 7 



+ 



9 2 T 

w 



+ 



3k 3T + 9k 9T + w 
9x 9x 9y 9y l 



9T 

P C 9t 



( 1 ) 



in the steady state system, the origin will be defined to 
be centered at the moving arc. Hence, the following trans 
formation-equations apply: 



x = x +v 
t = € 

where the primed quantities refer to the new system. 



Hence , 

9Jc = 9x , • 9Jc = -v9t , 9€ = 9t , = 0 



with these partial relationships, derivatives in the new 
system become, 



9T _ 9T %x + 9T 9T 

9x M 9x 9lf 9 x 9"x 



9 2 T _ 9 2 T 
dx* ? ~dx T 

9T dx 9T 9T 9€ _-v9T 9T 

9 t 9t 9t + 9t 91: 9x 9? 



Since temperature is constant in the steady state system 



9T 

9€ 



0 



and 



9T ~ V 9T 

9t dJ 
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Substituting into equation (1) and dropping prime notation, 



9 2 T 9 2 T 1 9k 9T 

9x2 9y2 k 9x 9x 



9k 9T 
9y 9y 




-v 9T 
a 9x 



where 



a = k/pc 



The equation is simplified by noticing that 



1 9k 9T = 1 9k ,9T? 
k 9x 9x k 9T . 9x' 



Finally, 

9 2 T 9 2 T 1 9k 

93F dy* k 9T 




-v 9T 
a 9x 



APPENDIX 2 
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DERIVATION OF HEAT LOSS TERMS 



A. Nucleate Boiling 



The Rohsenow 



15 



correlation 



is 



C L^ T w“ T SAT^ 



fg 




(g/A) N 

U L h fg 




1.7 



All quantities are taken at saturation (14.7 psia, 212°F) 
from the CRC HANDBOOK OF TABLES FOR APPLIED ENGINEERING 
SCIENCE . 

C T = 1.005 BTU/lb-°R 

Jj 

Pr L = 1.9 

u t ■ = 62.5 x lo " 7 lb-sec/ft 2 

JJ 

V 9 = 1 

hfg = 970.3 BTU/lb 

p T =59.8 lb/ft 3 

J_i 

p v = 0.0374 lb/ft 3 

a = 40.4 x 10 ^ lb/ft 

Substitution of these values into the nucleate boiling 
equation yields, 

% = U.542 (T w -T SAT) ] 1 , BTU/hr-ft 2 . 

in normalized form, 

q* = [1.542T* (T m -T sat )) 3 



where , 
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T* 




T 

SAT 

T 

SAT 



B. Film Boiling 



The Borishansky relation reported 



q F 



0.425 



K^gp^pL-Pv* (h v" h L } 
u v [V/g(p L” P v )]i 



by 

g 1 
3 



Jordan 



25 



is 



(T -T 
w 



SAT 




where , 

P L = 59.8 lb/ft 3 

p^< = 0.0374 lb/ft 3 

a = 40.4 x 10 -4 lb/ft 

h f = 970.3 BTU/lb 

u v = 2.5 x 10 7 lb-sec/ft 2 

K y = 14.3 x 10 _3 BTU/hr-f t-°R 



Substitution of these values into the film boiling 
equation yields, 

q F = 32.4 (T w -T sat ?* 75 , BTU/hr-f t 2 . 
where the normalized form becomes, 

q* = 32.4 [T*(T m - T sat )?' 75 



C. Transition Regime 

The proposed relation for heat flux in this regime is 

q = M(T -T ) k 
m T v w sat' 

to which the boundary conditions are applied: 

1) q T = q N at T w = 300°F 

2) q T = q F at T w = 500°F 



